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We investigate a cosmological model with a big-brake singularity in the future: while the first 
time derivative of the scale factor goes to zero, its second time derivative tends to minus infinity. 
Although we also discuss the classical version of the model in some detail, our main interest lies in 
its quantization. We formulate the Wheeler-DeWitt equation and derive solutions describing wave 
packets. We show that all such solutions vanish in the region of the classical singularity, a behaviour 
which we interpret as singularity avoidance. We then discuss the same situation in loop quantum 
cosmology. While this leads to a different factor ordering, the singularity is there avoided, too. 

PACS numbers: 04.60.Ds, 98.80.Qc 



I. INTRODUCTION 

It is a well-known fact that general relativity is an in- 
complete theory in the sense that solutions to Einstein's 
equations can contain singularities. These are regions 
(outside spacetime) where the theory itself breaks down. 
According to the singularity theorems, the occurrence of 
such singularities is a generic feature of 'physical' solu- 
tions to Einstein's equations. 

One outlook on this problem is to consider a quantum 
theory of gravity as the necessary completion of general 
relativity Consequently, it is expected that such a 
quantum theory of gravity is in some sense (still to be 
specified) free of these singularities. Investigations to 
this end are usually carried out, not in the full quan- 
tum gravity candidate theories, but in reduced models. 
That is, one takes a specific solution (conventionally and 
pragmatically specified by some symmetry) to Einstein's 
equations and in some way applies the quantization pro- 
cedure of the full theory to the reduced model. ^ Proto- 
types for such symmetry-reduced models are black-hole 
spacetimes and cosmological spacetimes. 

In our paper we restrict the discussion to cosmological 
models. Here, in the canonical approach, we have basi- 



^ A counter-example is causal dynamical triangulation [J, Q . Here 
exists the possibility to reduce the full quantum theory by inte- 
grating out all degrees of freedom except the scale factor. The 
resulting theory yields an action differing from the geometrody- 
namical minisuperspace action by an overall minus sign in the 
realm where the continuum limit is valid. Numerical evaluation 
predicts a closed universe undergoing a bounce upon reaching 
small scales. Moreover, quantum spacetime on these scales is 
predicted to be of fractal structure and dimension 2, coinciding 
with results obtained in the asymptotic-safety approach [l|, . 



cally two candidates for a quantum cosmological theory: 
minisuperspace quantization in the framework of the ge- 
ometrodynamical approach and loop quantum cosmol- 
ogy 0,11 i. In both approaches, one has to investigate 
whether singularities 'do not occur'. This implies that for 
each approach one has defined what the sentence 'singu- 
larities do not occur' means. To come to the point, for 
neither of the two theories a strict proof of the avoidance 
of singularities exists. 

Both approaches describe the universe via a wave func- 
tion on configuration space which has to be the solu- 
tion of a constraint equation. The constraint equation 
is the quantized version of the Hamiltonian constraint. 
The difference between both approaches lies in the way 
this equation is quantized. In loop quantum cosmology, 
one uses a so-called polymer representation instead of the 
conventional Schrodinger representation. This is done in 
analogy to the full theory. This procedure carried out 
in a naive way, leads to a difference equation in steps of 
a smallest length /ig- In geometrodynamics, one arrives 
at a differential equation, the Wheeler-DeWitt equation. 
In the continuum limit, /xq ^ (suitable conditions on 
the higher derivatives of the wave function implied), the 
loop quantum cosmological difference equation fades into 
the Wheeler-DeWitt equation [6|]. 

Recently, Ashtekar et al. [tI, 0| extended the ansatz 
using /io, replacing it by /I, which is a function of the 
densitized triad operator p. The equation is then a dif- 
ference equation in eigenvalues v of the volume operator, 
and the Wheeler-DeWitt equation follows in the con- 
tinuum limit for large volume. The factor-ordering of 
the Wheeler-DeWitt equation then does depend on the 
factor-ordering chosen for the difference equation. In Q 
and 0, S] different factor-orderings have been chosen. 
The two difference equations, in /io or /l, can be under- 
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stood in a broader context as implementing different ac- 
tions of the full Hamiltonian constraint. They are thus 
just two special cases of a wider class of constraints that 
might arise, the actual form of which should in principle 
be determined by the full Hamiltonian constraint, ^. 
Whereas in the first case, the coordinate edge length of a 
holonomy is fixed and does not depend on the scale fac- 
tor, in the second case it does. This can be interpreted 
as an implementation of the fact that in the full theory, 
the Hamiltonian constraint (whatever its exact form may 
be) creates vertices (in addition to changing the edge la- 
bels of the existing edges). As new vertices are created, 
the edge lengths decrease. The altered dynamics using 
p. then corresponds to a lattice in which the number of 
vertices grows linearly with volume. 

In loop quantum cosmology, results on singularity res- 
olution fall into one of three categories, As a first 



result one may quote that, in the isotropic case, the evo- 
lution equation is well-defined also on an evolution across 
the singularity. This is due to the discreteness of the evo- 
lution parameter which is a feature inherited from the 
full theory through the use of the polymer representa- 
tion, [13] ■ This allows to evolve a wave packet, starting 
from any initial state, deterministically across the singu- 
larity, [li]. 

A second hint on singularity avoidance, so far studied in 
isotropic models with massless scalar field 0, curvature 
index /C = 0, 1 and zero as well as non-zero cosmological 
constant, is the occurrence of a so-called 'bounce'. As a 
bounce one describes a deviation from the classical be- 
haviour such that a semi-classical wave packet starting 
on a classical trajectory for large scale-factor deviates 
from this trajectory upon approach of the classical sin- 
gularity and instead avoids the region of configuration 
space where the singularity is located. Here, avoidance 
refers to an exponential fall-off (in 0) of the wave func- 
tion, @,0,i. 

A third criterium is given by the boundedness of the ex- 
pectation value of the operator corresponding to the in- 
verse scale factor. As the inverse scale factor is related 
to the curvature in isotropic, homogeneous models, this 
hints at avoidance of the curvature singularity. This is a 
feature which follows from the use of holonomies as basic 
variables. It is a purely kinematical result as the expec- 
tation value is evaluated with respect to states from the 
kinematical Hilbert space, [HQ. 



The robustness of these results is disputable to differing 
degree. Whereas the possibility to evolve the wave packet 
through singularities in a well-defined way seems to per- 
sist in the full theory, this is not so clear for the other 
two criteria. 

The boundedness of the inverse scale factor seems to 
carry over to the full theory only when evaluated on a 



subspace of the kinematical Hilbert space, More- 
over, the occurrence of a bounce seems to be knit to 
isotropic models and even there it is not clear whether 
it should persist for more general settings involving a 
matter potential. The underlying concept in the models 
studied in this context is to use the scalar field as a 'time' 
variable (emergent time) with respect to which the wave 
packet is evolved (numerically). Transferring this con- 
cept to more general models including a scalar-field po- 
tential, one has to cope with a 'time' (i.e. (jj-) dependent 
evolution operator which is given by the square-root of 
the gravitational Hamiltonian plus the scalar field poten- 
tial energy. It can therefore be arbitrarily complicated. 
In addition to that, it is not clear that <j) defines a 'good 
clock' throughout the universe evolution. The advantage 
of this approach, on the other hand, is the existence of 
an inner product which is uniquely defined by a complete 
set of Dirac observables, and thus provides expectation 
values of observables. Most importantly, the inner prod- 
uct supplies the model with a probability interpretation 
(even though no connection to the measurement process 
is made). 

In the geometrodynamical framework, several mod- 
els have been investigated regarding their ability to re- 
solve the singularity problem. In this setting, singularity 
avoidance is defined as either a vanishing of the wave 
function at the point of the classical singularity^ or a 
spreading of semi-classical states denoting a break-down 
of semi-classical concepts in general (the end of the world 
as we know it). In the semi-classical regime, an approxi- 
mate Schrodinger equation can be derived from a WKB- 
expansion defining a notion of time [l[ . This time label is 
necessary to stack together the 3-hypersurfaces on which 
the wave function has support. The thus obtained 4- 
dimensional spacetime can now be probed for geodesic 
completeness. Only in semi-classical regimes a notion of 
geodesies exists, and thus we can speak of the existence 
of singularities — in the strict mathematical sense of the 
singularity theorems — only there. 

Accepting both criteria, singularity avoidance was 
found for big-bang/big-crunch singularities in various 
models (different scalar field potentials, cosmological con- 
stant, etc.) and for the big-rip singularity occuring at 
large scale factor, [15|. The big-rip singularity is a sin- 
gularity which the universe can encounter when it ex- 



^ More generally, it would be sufficient to demand that the proba- 
bility vanishes there; for example, the ground-state wave function 
for the hydrogen atom, as found as a solution to the Dirac equa- 
tion, diverges for r — > 0, but the probability there is zero because 
of the r^-contribution from the measure. In quantum cosmology, 
this question is more subtle because the fundamental measure is 
not known [l]]. 



3 



pands "too rapidly" 16|. This singularity occurs when 
the cosmological radius of the universe a{t) tends to in- 
finity at some finite moment of time simultaneously with 
its time derivative d{t) in such a way that the Hubble 
variable H(t) = a/a tends to infinity as well. Interest in 
this type of singularity is connected with the fact that it 
arises quite naturally in cosmological models with phan- 
tom dark energy, that is, dark energy such that the equa- 
tion of state parameter w = p/p < —1 17|, |18[, where p 
and p denote pressure and energy density of the cosmo- 
logical fluid, respectively. 

In the following, we want to analyze whether the so- 
called big-brake singularity can be avoided in a similar 
way. The big brake belongs to another class of cosmo- 
logical singularities not connected with the divergence of 
the Hubble variable itself but of one of its higher deriva- 
tives. Singularities of this type are called soft, quiescent, 
or sudden 19, [13, 21 1. These singularities occur at fi- 



nite value of the scale factor and its time derivative and 
hence of the Hubble parameter, while the first or higher 
derivatives of the Hubble parameter are divergent, which 
implies divergence of some curvature invariants. The big 
brake is a special example for a model from this class; it 
was first considered in [22] (see there the discussion after 
Eq. (2.13)) and later discussed in detail in [2l|. It can 
arise in tachyonic cosmological models 23| with a partic- 
ular potential: at some finite moment of the cosmological 
evolution the universe stops at finite value of its cosmo- 
logical radius with an infinite deceleration a — > — cxd. It 
was also noticed that the big-brake singularity can arise 
in more simple cosmological models, such as a universe 
filled with a perfect fluid obeying the equation of state 
p = A/ p, where A is a positive constant. This equation 
of state was considered in 24] in the context of wiggly 
strings (these are cosmic strings with small-scale wig- 
gles imposed on their dynamics). A fluid obeying this 
equation of state can be called "anti-Chaplygin" gas in 
analogy with the gas with Chaplygin equation of state 
p = —A/p, which has acquired some popularity in cos- 
mology as candidate for unifying dark energy and dark 
matter 2^ 26 1. Independent of the possible relevance of 
such a model for the real Universe, it has the merit of 
showing that intriguing features can occur in the quan- 
tum version, connected with the presence of a quantum 
phase at large (instead of small) scale factor. Quantum 
effects at large cosmological scales have previously been 
studied in the context of a classically recollapsing quan- 

ii. 
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tum universe 

Our paper is organized as follows: In Sec. II we present 
a simple classical model exhibiting a big-brake singular- 
ity. In Sec. Ill the Wheeler-DeWitt equation for this 
model is studied and approximate solutions describing 
wave packets are found. Their behaviour demonstrates 



that the classical singularity is avoided. Sec. IV contains 
a discussion of the big-bang singularity. Sec. V makes 
a comparison with the description of this model in loop 
quantum cosmology. Sec. VI contains a discussion and 
an outlook. Some technical details are relegated to an 
appendix. 



II. THE CLASSICAL BIG-BRAKE MODEL 

We consider a flat Friedmann-Lemaitre universe filled 
with a perfect fluid mimicked by a homogeneous scalar 
held. We require the fluid to obey an "anti-Chaplygin" 
equation of state p — A/ p, where p is the fluid pressure 
and p its energy density. In terms of the scalar field, 
these read 



The corresponding action is 



(1) 



S 




(2) 



where n"^ = 8ttG, N is the lapse function, A the cos- 
mological constant, V{(l)) a potential of the field 0, and 
K. = 0,±1 is the curvature index; we set c = 1. Fur- 
thermore, we set iV = 1, so the time parameter is the 
standard Friedmann cosmic time. The action then be- 
comes 



S ^ ^ f dt {-ad^ +JCa- ^a^) 



-y di (a302 _ 2a3y(0)) 
The canonical momenta are given by 



6aa 



(3) 



(4) 



The canonical Hamiltonian 7i, which is constrained to 
vanish, reads 



n tt; a 3/Ca 



0. (5) 



In the following, we restrict the analysis to fiat cosmolo- 
gies, /C = 0, without cosmological constant, A = 0. The 
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Hamiltonian constraint yields the Friedmann equation 
k2 k2 / ^ 



(6) 



3 3 \ 2 
The fluid obeys a continuity equation, 

p = -2.H{p + p) , (7) 
which in terms of the scalar field reads 

(8) 



• dV^ 
3H<j) + — = 



Using the equation of state, p = A/p, ([7]) can be easily 
solved for p in terms of the scale factor, 



p{a) 



B 



A 



(9) 



where i? > is some integration constant, and we have 
chosen the solution with p > 0. Note that p is well defined 
only for a < = (B/A)^^^, cf. Figure 1. As a* is 
approached, the density goes to zero. We note that B 
has dimension mass squared, and A has dimension mass 
squared over length to the sixth power. 
Using the result 0, one gets from ([S]): 



da 



Ah^) 



V3 



(to-t) , 



(10) 



where a{to) = a* ("big brake") and a(0) — ("big 
bang"). In order to calculate this integral, we substi- 
tute z = {B/a^ - Ay/'^, with < z < oo. Then (fTO)) 
becomes 



d^ . , — 
z^ + A 2 



{to - t) . 



(11) 



The integral on the left-hand side can be found in [2£ 
For (fTT|) one then gets 



In 



z2_^l/4^V2 + Al/2 



^1/4^^2 



+2 arctan ■ 



+ Tr9{z^ - A^/^) j (12) 
(13) 



We have added the Heaviside 6'-function in order to make 
the arctan- function continuous at the point z"^ = A^l"^ . 

For the total time that elapses from big bang to big 
brake one then gets 



2 \-~ z 
to = I dz 



K^/3 Jo + ^ VGuA^/^ 



(14) 



The solution for a{t) is shown in Figure 2. A simple 
approximate solution can be found in the vicinity of a*. 
To this end, we write a = — Aa, which simplifies the 
above integral to 



Aa 



dAa 



1 



a^,(6Aa) i 



(t-to) , (15) 



yielding 



Aait) = [C{to~t)]-^ . (16) 
So we find for the scale factor and its derivatives 



a{to) = a*, a{to) = 0, d{to) 



(17) 



At to, the evolution of the scale factor comes to a halt. Its 
'speed' is zero due to an infinite negative acceleration. It 
is this peculiar feature that gave the singularity its name, 
big-brake singularity. 

The first and second time derivatives of the scale fac- 
tor in terms of the scale factor itself are given by simple 
expressions. To this end, note that pO|l can be differen- 
tiated with respect to a, thus connecting d{t) with the 
scale factor according to 



da 



B 



(18) 



cf. Figure 3. Obviously, as a — s- a^,, a ^ 0. Differentiat- 
ing again with respect to time, one finds 



d^a fB 

d^-y"'^-^ 



1 - 



B 
4^ 



^-A 



(19) 



showing that d{t) — oo as a — > a^,, cf. Figure 4. 

What remains to be found, is an equation for (p. As 
we are interested in the quantum model, the solution in 
configuration space, (/"(a), suffices. This is obtained from 



= P + P , 



(20) 



using the equation of state and the Friedmann equation 
(l6|). The (exact) solution is 



^^(a) =T\/^artanh(^/l-^ ) , (21) 



B 



cf. Figure 5. This is only consistent if the potential is 
chosen to be 
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FIG. 1: Evolution of the en- 
ergy density p of the scalar 
field with scale factor a. 



FIG. 2: Evolution of the 
scale factor over cosmic 
Friedmann time t. 



V{(l>) = Vo I sinh(V3^|(/'|) - - 



sinlif VSk^I^I 



(22) 



Given the trajectories 4>{a) and a(t), the latter in explicit 
form only in the vicinity of the singularity, the classical 
model is thus fully described. Note that Vb = ^JA/A. 
From (fT4)) we find for the total lifetime of this model 
universe the expression 



io « 7 X 10" 



1 



(23) 



This lifetime is much bigger than the current age of our 
Universe if 




FIG. 3: Dependence of the 
derivative of the scale factor 
on the scale factor itself. 



Vo < 2.6 X 10" 



-30 



cm-^ 



which is, of course, a reasonable result because the crit- 
ical value of Vb just corresponds to the scale of the ob- 
served dark-energy density. 



III. THE QUANTUM BIG-BRAKE MODEL 

A. Wheeler— DeWitt equation 

Quantization is carried out in the canonical approach. 
Implementing the Hamiltonian constraint via Dirac's 
constraint quantization, one arrives at the Wheeler- 
DeWitt equation 




FIG. 4: Cosmic accelera- 
tion depicted over a. 
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FIG. 5: Classical trajectory 
in configuration space. 




= 0, (24) 

where a = Ina and the Laplace-Behrami factor order- 
ing has been employed. As we are interested in the be- 
haviour in the vicinity of the big-brake singularity, where 
4> is small, it is sufficient to approximate the potential 
there. We find 

2 \ 6 da^ d(t>^) ^ \<t>\ ^ ' 

(25) 

where Vq = Vo/3k^. 



normalizable solutions are given by 

J: ^ 1 

fk{xk) = NkXke~~l^k-i{^k) , (27) 

where Xk = '2.\J — and h\_^{xk) denote the as- 
sociated Laguerre polynomials; Nk ~ l/fcs is the nor- 
malization factor; fc S N 

The choice of the normalizable solution to ((26| is en- 
forced through the condition on the wave function im- 
posed for large \(f)\, cf. Sec. IIVDI The exact normalizable 
solution to ([26|) with the exact potential possesses a dis- 
crete spectrum; coincidence with the behaviour at small 
|(/)| is thus only achieved if the normalizable solution ([27]) 
is selected because the non-normalizable solutions have 
a continuous spectrum. 

Note that (pkixk) — * for 0, since L^_-^(0) = k. 

To simplify notation, introduce Z{a) = fi^/Va and Vq, = 
Voe^". Then, Xk = 2\(j)\/Z{a)k. The functions (pk{xk) 
are orthogonal such that^ 

J d4> Lpk{xk)<fi{xi) = Z{a)Ski ■ (28) 

The energy eigenvalues are 

i^.H = -j|.. (29) 

Inserting this ansatz in (j25p and carrying out a Born- 
Oppenheimer approximation, the resulting equation for 
Cfe(a) becomes 

^'^(^)-J^^k{a)=0, (30) 

where dots denote derivatives with respect to a. Thus 
Cfc is given by 



B. Born— Oppenheimer approximation to the 
Wheeler— DeWitt equation 

Equation (j25p can be solved, at least approximately, 
making the ansatz ^E* (a, 0) = C'fe(a)(^/c(Q:, 0), where 
(y9fc(a, 0) is the solution of 

- + ll^'") = Ek{o^)Mc^^^) . (26) 

cf. also [i^l, where a similar ansatz was made. We recog- 
nize that this is the radial part of the time-independent 
Schrodinger equation for a particle in a Coulomb poten- 
tial with I = and the wave function ripk- Thus, the 



where Iq, Kq denote modified Bessel functions of first and 
second kind, respectively. As a boundary condition, we 
require that the solution should vanish in the classically 
forbidden region, a > a^,. Therefore, ci — and only the 
MacDonald function Kg remains as solution. On the level 
of the Born-Oppenheimer approximation, the complete 
solution is therefore given by 



^ The validity of this relation is clear from the property of the ipn 
being eigenfunctions of a Hermitian operator; its direct verifica- 
tion is discussed in [sill . 
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k=l 



1 

Vc 



(2^I<^|) e-^LLi (^2^101 



(32) 



C. Derivation of classical equations of motion from 
the principle of constructive interference 

To derive a phase from this expression, approximate 
([26|) and ((30|) further by a WKB-approximation. Mak- 
ing the ansatz (^fe(a,(?!)) = e^'^^o^"''^) in Ck{a) = 
one obtains to zeroth order in h the 
Hamilton- J acobi equation for the (j)- and a-part, respec- 
tively. Integration yields for 5**0(0!, 0): 



arcsin 1 — 





TT 


h^^k^ ) 


" 2 




TT 


2ri2fc2 


~ 4 



(33) 



in which the Langer boundary condition at the a- 
dependent turning point (j)t {ot) = 2fi?k^ /Va has been em- 
ployed. From (f30|) . no phase results. This coincides with 
the limit h ^ in PT|) . as lim^j^oo Kq {x) « V^^"^- 
So ^^^(q!, 0) constitutes the entire phase. 

The classical equations of motion should follow from 
the phase through the principle of constructive interfer- 

0: 



9Sto 
dk 



k=k 



fcO I 



dk 



\k=k 



arcsin 1 



Vo.\<f>\ 



1 - 



2h^k^ 



, 



(34) 



Here, k 



Vo 



/ 3k, 



2 h 



This constant arises under the 



conditions that, firstly, k and so also k have to be dimen- 
sionless, and that, secondly, the only constants of the 
model are Vo (or Vo), (or A and B), h and k. With 
this choice, (l34l) simplifies to 



ds: 



dk 



^^k^k = h 



-arccos 1 




(35) 



For the classical trajectory, (pij) . this is 



dS' 



dk 



feO I - - ft 
k=k — " 



—arccos 1 — 



cosh 



cosh I 



2- 



cosh' 



(V3^|0|) 



(36) 



But the classical equation of motion was derived using 
the full potential. The quantum theory uses an approxi- 
mation to the original potential which is valid up to order 
O for small (j). Applying the same approximation 

to (15^. one finds 



dS' 



kO I 



dk 



\k=k 



(37) 



and so the classical solution pT|) satisfies the condition 
for constructive interference with the above choice for k 
for small (j), which is consistent with the approximation 
of the potential in 

There is, of course, also the question whether the 
Born-Oppenheimer approximation employed in the last 
subsection is a feasible approximation. We show in Ap- 
pendix A that this approximation is fulfilled in the limit 
a — > a*, which is just the region under consideration here. 



D. Singularity avoidance 

Wave packets in quantum cosmology have been con- 
structed in order to study aspects of the quantum-to- 
classical correspondence, in particular the validity of the 
semi-classical approximation [l[ [s^, 32 1. They are also 
useful in order to provide a consistent picture of the pre- 
big-bang to post-big-bang transition in quantum string 
cosmology [33|. such a construction is also useful in the 
study of singularity avoidance. 

Wave packets constructed from the solutions of ((25|) 
are of the general form 



00 

E 

k=l 



AkCk{a)'Pk{a, 4)) 



(38) 



We can choose initial conditions on a hypersurface a — 
ao- Here, it suffices to fix the values \l'(aoi</') s^nd 
^^da'^^ \ a=ao ■ fo^' ^hc choscu normalizablc solution 
([27)1 (fkia, (f) vanishes at </> = for all k and a, the wave 
packet is zero there. This is, of course, independent of 
the initial conditions. But the classical singularity occurs 
at = 0. So out of these solutions, no wave packet can 
be constructed which does not vanish at the classical sin- 
gularity. Taking a as an internal time variable, one can 
calculate the probability distribution, 
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I^P(ao» = '^AkAiCk{ao)Ci{ao)ipi{ao,(l))'fk{ao,(l}) , 

l,k 

(39) 

for each 'instant of time' ao. It is obvious that 
|4'P(q!o,0) = at = 0. We emphasize that this is 
a consequence of the choice of (p7)) . 

To manifest the ehmination of the classical singularity 
on the quantum level, also expectation values have been 
employed, see, for example, [7|. Before calculating the 
expectation value for for this model using the inner 
product ([28l) . recall that the avoidance of the singularity 
of the Coulomb potential in ordinary quantum mechanics 
is caused by a lowest bound on the energy due to quan- 
tization. This again leads to a minimal radius for the 
'trajectory' of the electron. 

Analogously to the Coulomb potential in ordinary 
quantum mechanics, the energy (of the matter compo- 
nent) in our model is also bounded from below. The 
minimal energy, given by (|29p for k = 1, corresponds to 
a minimal 'radius', that is, to a minimal value for |0|. 
This is given by 



We want ^{uq, (f) to be a Gaussian centered at 0o with 
width \J~^, where Zq = Z^ao). The center should be 
the value of the classical trajectory at ao. Note that we 
have two classical solutions, 0+ and 0_, see (PT|) ."^ So in 
fact, we have to construct two Gaussians, one centered at 
00, the other at — (/)o and superpose both. Write therefore 

*(ao, cj)) = *_(ao, + ci^'+(ao, , (41) 

where ^+ denotes the part of the wave packet being cen- 
tered around 0o and '^^ the part centered around ~(j)Q 
at initial 'time' ao. 

The calculation of the wave packet will employ only the 
WKB solution of ([30|) . With suitable initial conditions, 
it reads 

C,(a)- ('^Vexp[-i-^,/A(e«"-e«-)' 

(42) 

Introducing t = e^" (and denoting To = e^"°). 



(|0,|)(a) = [Cfe(a)]^^[Z(a)]^fc2 



f 1 Vc. 



3^ 

2^ 



for fc = 1. The classical singularity lies at a = a*. In 
this case the minimal energy is given by 



y2 



(40) 



and the expectation value for 101 is consequently given by 
(1 01 1) (a*). The boundedness of the energy here prevents 
the scalar field to evolve to the singularity, |0| = 0, in 
this case as well. 

Note that for a — > cx), the energy is no longer bounded. 
In this case (|0i|) — > 0, cf. (|40|) . Of course, one should 
keep in mind that the expectation value in quantum cos- 
mology has no interpretation in terms of measurement 
results as it has in conventional quantum theory. 



E. Construction of wave packets 

Apart from the avoidance of the singularity, we want to 
study semi-classical and quantum regimes of the model. 
To this end, we construct semi-classical wave packets and 
study their behaviour. Especially we are interested in the 
regions of configuration space where these packets spread 
(if they spread at all). 



exp 



1 Vo 
6 V2?i2fc2 V K 



(r-ro) 



(43) 



Start with the ^f^-part of the wave packet. We here find 
the requirement 



*+(ao,0) -^^+V'fe( 



= e ^0 . (44) 



k=l 



Decomposing the Gaussian into the tfk{oiQ,<t>), one ob- 
tains for the coefficients the somewhat lengthy expression 



^exp 



— i— 

2Zo \2k 



E™=o(-ir(^" + i^ 



(k\\ 



\\l ZOk) 



loi] D-(™+2) 



(fc-m-l)!(T?i+l)! 
hi^k-^^) 



(45) 



where I?m(x) denote parabolic cylinder functions. Note 
that this expansion in Lp^ cannot be performed at = 0. 
Here, <^k{oi, = 0) = for all k as remarked above. 

The amplitude of is obtained in a similar way (or 
by just substituting — 0o for 0o). The solution is 



^ The case with two Gaussians is the most general case. One may, 
of course, wish to choose only one Gaussian in order to represent 
only one branch of the classical solutions by a wave packet. 
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A: 



exp 



Zo 



^+2^(21 + ^0)' 



E™io(-ir(™ + i)Ti^ 



l)!(m+l)! 



/ / \ 

So the wave packet is given by 



(^ + 0o) . (46) 



(47) 



The total probabihty for the wave packet is calculated 
via 



exp 



(-i 



Vb 



(r-ro) 



(48) 



Probability is thus not conserved with respect to internal 
'time' T, as expected [1]. Choose the normalization of the 
wave packet such that at ao, / d0 j^tp = 1. Then, 



*(i-,< 




FIG. 6: This plot shows the wave packet. It follows classical 
trajectories with initial values ao = 1 and 00 ~ 0.88. The 
classical trajectories are depicted in the (r, 0)-plane; recall 
r = a«. 

This corresponds to a singularity occuring at a^, = 10 e. Pa- 
rameter values are Vo = 1, ft = 1 and Ci = 1, cf. (|49|l . 
Summation was carried out up to fc = 50. 



^ 00 



(49) 



fc=i 



where the normalization factor is given by 



one has \(j)\ — > 00. Thus one can approximate the poten- 
tial by an exponential in the vicinity of this singularity. 
Choosing units such that = 6, one obtains the follow- 
ing form of the Wheeler-DeWitt equation: 



C = 



(50) 



A plot of the wave packet is shown in Figure 6. We 
recognize that the wave function is peaked around the 
two branches of the classical trajectory in configuration 
space, but goes to zero if the region of the classical big- 
brake singularity, a — > a^,, is approached. In this sense 
the classical singularity is avoided in the quantum theory. 
This is a consequence of the choice of the normalizable 
solution ([27|) . which vanishes at = (the region of the 
big-brake singularity). Moreover, we find that the wave 
packet does not spread along the classical trajectory. 

IV. REMARKS ON BIG-BANG SINGULARITY 



Introducing coordinates zi — a + — a — |0|, the 

Wheeler-DeWitt equation becomes 



92 



dz\dz2 



^ = /(Z1,Z2)* 



(52) 



One can now find coordinates so that the function on the 
right-hand side cancels. One is left with 



92 



52 



where 



+ * = 



(53) 



A. Solution to the Wheeler-DeWitt equation 

So far, only the big-brake singularity of the model was 
considered. But the model possesses a second singularity. 
Namely, its evolution starts with a big bang: as a — > 0, 



w(a,(/)) 



2V% 3fa + ^|0|) 



coshX — sinhX 

sinhX ^ coshX 

V2 



,(54) 
,(55) 
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and X 



72^ 



A solution to this equation can 



D. Singularity avoidance 



be found from the WKB-ansatz ^' = / <lkA{k)e^Ti^°^ . 
Inserting this ansatz into ([53)1 yields the Hamilton-Jacobi 
equation of which an exact solution is given by S'ofc = 
ku-y/k^ - Iv. Of course, the Hamilton-Jacobi equation 
is also solved by actions with different signs in front of 
u and V. These are obtained from the one chosen above 
through rotations in the (u, w)-plane. As m > 0, only two 
solutions can be mapped onto each other. 

B. Recovery of classical trajectories 

The classical trajectory in the vicinity of the big bang is 
recovered using the principle of constructive interference 
i^^lj^^^ = 0. For A: = V2 one finds (j){a) = ±^a. This 
is just the classical trajectory obtained from (PT|) in the 
limit ^ 1 with initial condition B ~ ^ and fixed A. 

C. Construction of wave packets 

We get the following exact wave-packet solution to the 
Wheeler-DeWitt equation: 



^{u,v) ^ J dkA{k) (Cie-^('^-"-V(*='-i)'') +C2 c.c.) , 

(56) 

where c.c. denotes the complex conjugate of the precee- 
ding term. By construction, the classical trajectories can 
be recovered from this equation through the principle 
of constructive interference. Choosing as amplitude a 
Gaussian with width a centered around fc, 



A{k)^ 



and taking Ci = C2 for definiteness, one obtains wave 
packets of the form 



exp 



2ah 



1 - ia^hS^ 



2(cr- 



ihS" 



(57) 



where a Taylor expansion of S'ofc has been carried out 
around k (primes denoting derivatives with respect to 
k) and the terms of the order (fc — fc)'^ in the expo- 
nent have been neglected. (For simplicity, in this ex- 
pression S'ofc (fc) = So.) This can be done if the Gaussian 
is strongly peaked around fc, that is, if a is sufficiently 
small. Since SQj.(fc) — gives the classical trajectory, the 
packet is peaked around it. 



Due to the fact that u > 0, two inequivalent actions 
exist. Apart from the wave packet constructed from 
S'ofc = ku — ~ Iw, one gets a second wave packet 

-ku - 



ku — ~ 
constructed from Sofc = 



\Jk^ — \v. Moreover, the 
entire (a, 0) plane is mapped into only a quarter of the 
(m, v) plane. One would therefore require the wave packet 
to vanish on the boundary of the physical region. The 
only solution satisfying this requirement is naturally the 
trivial one. To get a non-trivial solution, one has to lessen 
the boundary condition and require ^' = only at the 
origin of the (u, v) plane. The fact that the wave packet 
does not vanish at the m = and v — Q line is due to the 
non-normalizability of the wave packet in both a and 0, 
which in turn has its origin in the fact that the approx- 
imation to the classical trajectory for large \(\)\ has no 
turning point. 

The implementation of the condition of normalizability 
results in a wave packet which vanishes at the big-bang 
singularity, >OasQ:— s-— 00, and spreads for large a. 
This is equivalent to the condition '5 ^ as |0| ^ 00. 
The condition implied for large |0| in the vicinity of the 
big bang thus implies and justifies the normalization con- 
dition imposed in the derivation of the solution to the 
Wheeler-DeWitt equation in the vicinity of the big brake, 
cf. (pS)) . We thus impose basically two conditions on 
the wave function. The first one is that — > when 
|(/)| 00, resulting in a normalization condition for the 
approximate solution in the vicinity of the big-brake sin- 
gularity and the elimination of the big-bang singularity. 
The second condition is to require ^E" — > as a — > cxd to 
ensure the existence of wave packets that follow the clas- 
sical trajectory. Upon matching the wave function in the 
two regimes, one would expect quantization conditions 
as observed e.g. in [s^l or [s^l. The big-bang singularity 
does therefore not exist in the quantum theory. 
The method employed in this section mirrors the calcu- 
lation carried out in [T^ . The picture one obtains is thus 
the following. For large |0|, the wave packet vanishes and 
so does the wave packet for small \(^\. In the intermedi- 
ate region, the packet is peaked around the corresponding 
approximation to the classical trajectory, cf. Figure 7. 



V. RELATION TO LOOP QUANTUM 
COSMOLOGY 

As discussed in the introduction, there is a whole class 
of difference operators in the current loop quantum cos- 
mology literature. The ambiguity stems not only from 
the freedom to choose a factor ordering, but also from the 
fact that the Hamiltonian constraint contains a curvature 
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FIG. 7: The difFerent regions of the wave packets and the 
classical trajectories are shown. 



term which, when expressed in terms of holonomies, is 
given by a hmiting procedure. This limiting procedure 
consists of shrinking an area to zero. But in loop quan- 
tum gravity there is a smallest area and thus the limit is 
heuristically reduced to 



Area ^ A = minimal area 



(58) 



There are at least two ways to implement this. First, 
one can send the side length of the area to the value 
/i ^ /io = [^. Then /io is just some number. But on 
the other hand, one can require that the physical length 
is taken to its minimal value. But the physical length de- 
pends on the scale factor and so does /i, — > /i = /i(|p|), 
\p\ = 0-^ Depending on which of the two viewpoints 
is taken, one arrives at a difference equation either in 
eigenvalues of the triad, /i, or in eigenvalues of the vol- 
ume operator, v. (As discussed in the introduction, this 
is more suitably understood as a volume-dependent cre- 
ation of vertices an thus a refined implementation of the 
action of the full Hamiltonian constraint.) 
The Wheeler-DeWitt equation is recovered in the respec- 
tive continuum limit. Starting from the same factor- 
ordering of the difference equation, both versions of it, 



the one equidistant in /x, the other equidistant in w, have 
the same Wheeler-DeWitt limit, meaning they yield the 
same factor-ordering of the Wheeler-DeWitt equation. 
As in and 0, [1] different factor orderings have been 
employed, we comment briefly on both of them here. 
The question of whether the preceding result persists in 
loop quantum cosmology can then be reformulated as 
the question whether the results obtained in Sec. IIIII are 
robust with respect to a change of factor ordering. 



A. Non-covariant factor ordering 

The question here is under which conditions the con- 
tinuum limit is justified. It is justified if the discreteness 
of spacetime is negligible compared to the length scales 
occuring in the model. For large scale factor, a 3> /io, 
one can argue that the limit /ig ^ is a sensible approx- 
imation. Thus singularity avoidance for large-scale sin- 
gularities as, for example, the big rip or big brake, in the 
loop quantum cosmology framework, reduces to singular- 
ity avoidance induced by the Wheeler-DeWitt equation. 

The Wheeler-DeWitt equation emerging in the con- 
tinuum limit of the difference equation employed in 
is 



fi2 
2 



k2 2 9^* a^^-i 

6 da?- d(p 



a^-^vl/ ^ , (59) 

m 



which differs from (pS)) by the choice of factor-ordering. 
Making the ansatz ^(a, 0) = A(fc)Cfe(a)(pfe(a, 0) and 
requiring iy9fe(a, (p) to be a solution of 



f^2 q2 y 

y ' = -~Ek{a)ipk{a,(j)) , (60) 



one finds as before the solution 



(61) 



% 2 



where Xk = 2y and Ekia) = -^^-^p-^o 

Then the equation for Ck{a) is given by 



d^Ckja) 
da2 

which is solved by 



a'"Ck{a)^0, 



(62) 



^ The interpretation of the area operator in loop quantum gravity 
is still unclear. The operator itself is not a Dirac observable 
and thus even in the sense defined by the loop quantum gravity 
community not an observable (though it can become a Dirac 
observable when matter is added). In how far the area operator 
relates to a physical area is thus unsettled. 



Ck{a)= c,^J^iy-^a' 



(63) 
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The complete solution has an analogous form to the 
quantum geometrodynamical formulation in Sec. Ill B. 
The decisive result is that, because only the factor order- 
ing of the gravitational part has been changed compared 
to P5|) . the solution for (^^(0, a) handles the singularity 
avoidance in this framework as well. 



B. Covariant factor ordering 

The factor ordering in the more recent paper f?., 
yields the Laplace-Bcltrami factor ordering for the 
Wheeler-DeWitt equation in the continuum limit. As 
this is the factor ordering we employed throughout this 
paper, the results of the previous sections carry over to 
the loop quantum cosmology analysis without alteration. 
Note, though, that a consistent loop quantization re- 
quires a polymer representation of the matter fields as 
well. This would require a Bohr compactification of 
which may bound the approximate potential V{(p) = 
from above. As the vanishing of the wave function at 
= is related to the divergence of the potential at this 
point, it is not clear whether the previous results would 
survive in the polymer representation; namely, it is imag- 
inable that the regularity condition and thus the ensuing 
condition that ipk{4> = 0, a) = becomes redundant. 
This has to be investigated in future publications. 



mology. Here, the analysis was restrained to the vicin- 
ity of the big-brake singularity. Two different factor- 
orderings were studied. For both we could corroborate 
avoidance of the big-brake singularity. 
Due to the special form of the potential, we were able 
to solve the model in the geometrodynamical framework 
also in the vicinity of the big-bang singularity. The choice 
of variables enforces a boundary condition which causes 
the wave function to vanish at the big bang. This sin- 
gularity is thus also eliminated in the quantum theory. 
The imposition of boundary conditions on both ends of 
the evolution, near the big bang and near the big brake, 
should imply some kind of quantization rule upon match- 
ing the wave packets in both regimes. Such a matching 
has not been carried out. 

What are the implications of this singularity removal? 
Since the wave packet starts to spread when approaching 
the region where the classical singularity would lurk, this 
means that the end of the classical evolution is reached. 
Any information gathering and utilizing system would 
stop to exist. A similar scenario may happen when the 
turning point of a classically recoUapsing quantum uni- 
verse is approached [27| . Classical time then comes to an 
end. The details of such a scenario can, of course, only be 
discussed if one goes beyond minisuperspace: the treat- 
ment of concepts such as entropy and the arrow of time 
need additional degrees of freedom [l|, . We plan to 
return to this issue in a future publication. 



VI. DISCUSSION AND OUTLOOK 

We studied a Friedmann-Lemaitre model with a scalar 
field obeying an 'anti-Chaplygin' equation of state. This 
model classically ends with a big-brake singularity. The 
singularity stands out because of its negatively diverging 
second derivative of the scale factor. This works as an in- 
finitely strong 'brake', forcing the derivative of the scale 
factor to go to zero. The evolution of the scale factor 
stops. Upon quantizing this model in the quantum ge- 
ometrodynamical framework, we are led to the Wheeler- 
DeWitt equation. It can be solved in the vicinity of 
the big-brake singularity. A separation ansatz yields a 
Schrodinger-type equation for the hydrogen atom for 
(which here plays the role of the radius in the quantum 
mechanical equation). Solutions to this equation vanish 
at (/) = 0, which corresponds to the singularity. Thus, 
independent of the choice of initial conditions, whatever 
wave packet is constructed out of these solutions, it is 
condemned to vanish at the singularity. Therefore we 
can conclude that in this model as well, the large scale, 
soft, future singularity is removed from the quantum the- 
ory. 

The same model was also studied in loop quantum cos- 
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APPENDIX A: VALIDITY OF 
BORN-OPPENHEIMER APPROXIMATION 

The Born-Oppenheimer approximation consists in ne- 
glecting cross-terms of the form 

f°° d 
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To be able to give some indication of the quality of the 
approximation, it is necessary to evaluate these terms. 
Carrying out the differentiation, one finds 



This integral can be evaluated using the general formula 
34| 



dxe-'"'x''L''JXx)L'^{px) 



r (to + n + a + 1) (6 - A)" (6 - /i)'" 

m\n\ ^rn+n+a+l 

Starting with the first part of the integral in (jAl[) 



d0 cf'e-'-*Li^^{\4>)L]M) 



where the short hands k = , A — . and a 

Z{a)n7n ' Z{a)n r~ 

— have been used, one finds 

Z{a.)m ' 



Xi = (-1)""^ ^i^m'^n{m + l)^{n + ni - 1)!(to + 1)! 
Here, 



Xn 

Xl2 
Xi3 



^ , ni — n \ I I m + n 

2 I I F I —TO, —n + 1; —to — n\ 

ni ~ n 

2\ 



_F ^— m + 1,— 7^+l;— TO — n + 1; 

(n + TO + 1) (to — n)^ 
to(to + 1) (to + n) 

F —TO — 1, — n + 1; —TO — rt — 1; 



m + n 



m — n 



m + n 



m — n 



Similarly, the second part. 



X. 



can be integrated to 



X2 = {-ir 



Z{af 



i{m - 1)2(to - 2)(to + 71 - 3)! 









^ riTO ^ ^ 


( TO 


— n\ 


\m — n 1 


V TO 


+ n / 



rn+n 



X [-^^21 — -^22 — X; 



23J 



(to + n— 2) f m — n 



where 



(to — 2) \TO + 

^— TO + 2, — n + 1; —TO — n + 2; 

X22 = —TO + 3, — n + 1; —TO — 71 + 3; 



m + n 



m + n 



X23 — 



(771 + 71 — 2) (to + 77 — 1) f m — n 



{m — 2) (to — 1) \ TO + 77 

—TO + 1, —77, + 1; —TO — 77 + 1; 



TO + 7T, 



m — n 



Taking into account the prefactors, one finds 



Anm = ^(a)3TO 



Xi 



x. 



So Anrn Z (a) cx fl^^. For a = a^, it thus takes its 
minimal value. This shows that the Born-Oppenheimer 
approximation is fulfilled best when one approaches the 
region of the classical singularity. 
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